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A new test of a wide class of interest rate models is proposed and applied to a recently developed 
quantum field theoretic model and the industry standard Heath-Jarrow-Morton model. This test is 
independent of the volatility function unlike other tests previously proposed in the literature. It is 
found that the HJM model is inconsistent with the data while the quantum field theoretic model 
is in significant agreement with data. We also show that a portion of the spread between long and 
short term interest rates is explicable in terms of this model. 
Physicists have been working on several aspects of financial research over the last decade Q. One of the most 
important and as yet unsolved problems in finance is the modelling of interest rates. The interest rates at any point 
in time form a usually continuous curve (current interest rates for different times in the future) called the forward rate 
curve (FRC). We denote these rates by f(t,x) where t is the current time and x is the time in the future for which 
the forward rate applies. For example, /(l, 2) is the interest rate one year from now for an instantaneous deposit to 
be made 2 years into the future. 

The earliest interest rate models (eg., Vasicek dealt only with the spot rate (current interest rate for the present 
time) and the forward rate curve was treated as a derived quantity. These were found to be inconsistent with the 
observed data. The Ho-Lee ||| and HJM Q models were developed to deal with this problem by modelling the entire 
FRC rather than just the spot rate. The HJM model is, however, limited in the sense that the Brownian motions on 
which the HJM economy depends are independent of x. One way of removing this restriction is by formulating the 
forward rate curve as a quantum mechanical string as in Baaquie Q (we will refer to this model as the quantum field 
theoretic model). 

Several empirical tests of the HJM model have been performed (eg., Biihler, Uhrig-Homburg, Walter and Weber 
||, Flesker J?J, Sim and Thurston g) with mixed results. All of the tests assume a certain form for the volatility 
function a. In this paper, we propose a test which is independent of the volatility function. The test is applied to the 
HJM and quantum field theoretic model which was introduced in j5j. The HJM model can be formulated as a limit 
of the model in |5j which is briefly reviewed below. 

The one factor quantum field theoretic model models the forward rates as 

df(t, x) = + 
ot 

where A(t, x) is a quantum field whose action is given by 

/•oo pt+T FR 

S[A] = dt dxC[A] (2) 

J t Jt 



(3) 



where Tfr is the largest time to maturity for which the forward rates are defined (the domain is hence a semi-infinite 
parallelogram defined by t > to, t < x < t + Tptf). Tfr is introduced to ensure that the action is well defined but does 
not affect final results as the limit Tfr — -> oo must be taken. When fi — > 0, this model reduces to the HJM model 
upto a rescaling (for details, please see || ) . 

We assume that the function a(t, x) depends only upon the variable 9 — x — t. This is a theoretically reasonable 
assumption as it is the result of assuming that the theory is time translation invariant. Most of the functions used 
for a(t, x) in the literature satisfy this condition. 

The initial forward rate curve f(to,x) has to be specified. The field values of A(t,x) on the rest of the boundary 
points of the domain are arbitrary and are integration variables. The presence of the second term in the action given 
in (J2j) seems to be justified from the phenomenology of the forward rates |^] and is not ruled out by no arbitrage . 

The moment generating functional for the quantum field theory is given by the Feynman path integral as 

If , /*'* dt f t+TFR dxJ(t,x)A(t,x) a, A] 

Z[J] = - DAe Jt o J* e s[A] (4) 

On performing the calculation (details are provided in [pj ) , we obtain 

i f* dt r FR ded9'.j{t,8)D{e,e';t,T FR )j{t,e') 
Z[J\ = e Jt o J0 (5) 
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where 9 = x — t, 9' = x' — t and the propagator D(9, 9' '; t, T) is given by 



D(9, 9'; t, T FR ) = ■ ■ 3 F ^ — - sinh /z(T FjR - 0) sinh /x0' 



sinh 3 (pT F R ) 

{1 + sinh 2 (^T Ffl )e(0 - 6>')} + sinh fi(T FR - 9') sinh pB 
{1 + sinh 2 (^T Ffl )e(0' - 9)} + cosh(nT FR ) 

{sinh /i0 sinh ji9' + sinh p(T FR — 9) sinh p(T FR — 0')} (6) 

The above results are for unconstrained boundary conditions. It is, however, well known that short term interest 
rates are heavily influenced by central banks. Hence, it is reasonable to treat the field at the boundary where t = x 
differently. If we assume that the field at that boundary (i.e. A(t, t)) is distributed normally with variance a, we 
obtain the propagator 

n (9 9>\-n(ft 9>\ MM m 

D 1 (9,9)-D(6,d)- D(m + a (7) 

(It can be readily seen that the mean of the field at the boundary does not affect any of the results due to the no 
arbitrage condition.) 

To understand the significance of the propagator D(9, 9'; t, T FR ) we note that the correlator of the field A(t, 9), is 
given by 

E(A(t, 9)A(t', 9')) = S(t - t')D(9, 9';t, T FR ) (8) 
It can be readily shown that the no arbitrage condition is satisfied only when 



a(t,x) = <r(t,x) J dx'D(x,x';t,T FR )a(t,x') (9) 
In the limit [i — > 0, D — > I and we obtain the well known result 

a(t, x) = a(t, x) f dx'a{t,x') (10) 



for the one factor HJM model. 

Following Bouchaud ||, we use the daily closing prices for eurodollar futures prices as a measure of the forward 
rates. The eurodollar futures prices are linearly interpolated to calculate the forward rates at 3 month intervals. 
The 3 month deposit rate that the eurodollar futures actually represents is taken to be a good approximation to the 
instantaneous forward rate. The data used for this paper are the same as that used in ||. The data cover the 1990s 
and the length of the dataset is 846 trading days and forward rates 7 years into the future are available. 

We parametrize the forward rates as f(t, 9) rather than f(t, x) as this considerably simplifies the analysis consid- 
erably since the domain shape in the (t, 9) variables is rectangular. 

We concentrate mainly on the following quantities (again partially following ||) 



v(9) - y/<sp(t,e)> (if) 

n(a , <Sf(t,9 mm )(Sf(t,9)-Sf(t,9 mm )) > 

C{9) <SP(t,9 mm )> (12) 

with the differences being taken over one trading day (e), 6f(t, 9) — f(t + e, 9) — f(t, 9) and 9 m i n being three months. 
We assume that there are 250 trading days in a year. In the following analysis, we use the discretization <5(0) = -. 

Using the one factor HJM model, we can derive the following expressions for the above quantities which are accurate 
to zeroth order in e 

V HJM (9)=a(9)V^ (14) 
C hjm {9) = JP--1 (15) 
tbjm{9) = (r(9 min )y/e (16) 



In deriving this equation, we have discretized the Brownian motion process W as W(t) = y ~x where x is a random 
number with the standard normal distribution. We particularly note that the ratio thjm{9) is independent of cr(9) and 
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is in fact constant. The ratio as calculated from the data is shown in figure I and can be seen to be far from constant. 
Hence we see that the time translation invariant one factor HJM model is inconsistent with the real evolution of the 
FRC for any choice of function <r(9). 

Using the unconstrained quantum field theoretic model, we can again derive the expressions for the above quantities 
to zeroth order accuracy in e to obtain 



V QFT (0) = <r(eWD(0,e;t,T FR )e 

a(6)D(6,6 min ;t,T FR ) 

<^QFT{V) 



(17) 
(18) 



Observed and Fitted Ratios 
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FIG. 1. The observed and fitted r(9) for the quantum field theoretic and HJM models. It can be seen that the quantum 
field theoretic model is in much better agreement with the data. 





/"(year x ) 


CT(6» TOm )(year _1 ) 


Least Squares 


0.0822 


0.0308 


Bootstrap (90% C.I.) 


(0.080, 0.085) 


(0.0299, 0.0317) 


Serial Segments (90% C.I.) 


(0.081, 0.099) 


(0.028, 0.033) 



TABLE I. Results obtained for the unconstrained quantum field theoretic model. 
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The ratio r(9) is thus given in this model by 



t qft(9) 



cr{Omin)\JeD{0, 9; t, TFR)D{6 min , 9 m i n ; t, T 1 



FR) 



D( 



\ 9 m in', t, TpR 



(19) 



which is still independent of cr(9). However, we note that the ratio is no longer constant. Using this fact, we can fit 
the ratio to find /x and <j(9 m i n ). We took the limit Tfr — > oo as required and used the Levenberg-Marquardt method 
p| to obtain the non- linear least squares fit. The results are shown in table |. The confidence intervals were obtained 
through the bootstrap method ]l2|] . An alternative confidence interval was obtained by dividing the data into series 
of 500 days starting from the first day, second day and so on and so forth. The function r(9) was calculated and the 
parameters fitted for the resulting 346 data sets. The confidence interval using these serial data sets are also shown 
in table Q. The fitted ratio together with the observed values are shown in figure 

Sigma from C and V 

for the QFT model 




FIG. 2. <r((9) derived from V and from C for the quantum field theoretic model. 



Sigma from C and V 

for HJM 




FIG. 3. <j(6) derived from V and from C for the one factor HJM model 



We can also use equations |l7] and |Ts| to obtain two different estimates of the function a{9). These two estimates are 
plotted in figure 0. Similar estimates of a(9) for the one factor HJM model are plotted in figure |. As can be readily 
seen, the HJM model is seen to be inconsistent with the data. On the other hand, the quantum field theoretic model 
is in good agreement with data. It is also interesting to note that the volatility function for the HJM model derived 
from the data is very far from the constant or exponential forms that are commonly used in the literature. 

Performing the same procedure for the constrained quantum field theoretic model, we obtain the results in table 
p. The fitted ratio in this case is shown in figure || The two estimates of a{9) are shown in figure |. The agreement 
between the two functions is better than in the case of the unconstrained model as may be expected due to the 
additional parameter involved. However, it can be seen from the large confidence intervals that the model is probably 
overspecified since different values of the parameters give rise to very similar values for r(ff). 

Another quantity that is of great interest is the mean spread between the forward rates and the spot rate 



s{6) =<f(t,9)-f(t,9 mm )> 



(20) 
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The spread is the linear sum of two parts : the spread due to the market price of risk and the spread that results 
from the no arbitrage condition in the model. Since we assume that a is only a function of 9, it follows that a is also 
only a function of 9. To calculate the spread due to the model, we assume that the initial forward rate curve is flat 
or that the effect of the initial forward rate curve becomes negligible after a long time. 



Observed and Fitted Ratios 

Constrained Boundary Conditions 




FIG. 4. The fit for r(6) obtained using the constrained quantum field theoretic model 



In that case, the mean spread due to the no arbitrage condition in the quantum field theoretic model is given by 

sqft{9) = (9- 9 min ) lim a(t) - ( a{t)dt (21) 
t_> °° Je min 

where 

a QFT (t) = a(t) [ a{6)D{t,6;t,T FR )d6 (22) 
Jo 





Myear x ) 


^(^nXyear- 1 ) 


a (year 1 ) 


Least Squares 


0.0174 


0.181 


0.0024 


Bootstrap (90% C.I.) 


(0.006, 0.020) 


(0.158, 0.548) 


(0.0002, 0.0032) 


Serial Segments (90% C.I.) 


(0.019, 0.055) 


(0.044, 0.203) 


(0.002, 0.071) 



TABLE II. Results obtained for the constrained quantum field theoretic model 
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FIG. 5. <y(9) derived from V and from C for the constained quantum field theoretic model 
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FIG. 6. The calculated spread due to no arbitrageand the observed mean spreads. The large difference is due to the spread 
due to risk aversion in the market. 

Using one of the estimates of <r(0) (using either gives very similar results), we can calculate the spread due to the 
no arbitrage condition by numerical integration. Due to the relative inaccuracy of the estimation of cr(6) in the first 
place, a trapezoidal integration was considered sufficient. The result together with the observed spread is shown in 
figure [| It is seen that the calculated spread is significantly smaller than the actual spread which is consistent with 
the existence of the spread due to risk aversion. However, we see that a significant portion of the spread might be 
derived from the way the forward rate curve evolves. A very similar result is obtained when the constrained quantum 
field theoretic model is used. 

To summarize, we have proposed a new way to test the one factor, time translation invariant Heath-Jarrow-Morton 
and Baaquie's one factor, time translation invariant quantum field theoretic model for the evolution of forward rates 
using historical eurodollar futures data. We have found that the one factor HJM model can be rejected while the 
quantum field theoretic model is consistent with the data. We also find that a quantum field theoretic model with 
constrained boundary conditions to reflect the special nature of the spot rate is also consistent with the data but the 
parameters of the model cannot be sufficiently accurately derived using this method. We also show that a significant 
portion of the spread can be explained by the quantum field theoretic model. 

We would like to thank Jean Phillipe Bouchaud for interesting discussions and Science and Finance for kindly 
providing us with the data. 
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